A COMPARATIVE STUDY OF STRESS - FIELDS OF BLADE - DISC 
ATTACHMENT CONFIGURATIONS OF A TURBINE BLADE 


PITAMBAR CHAUDHARY 


C 3930 



DEPARTMENT OF MECHANICAL ENOINEERING 

INDIAN INSTITUTE OF TECHNOLOGY, KANPUR 

OCTOBER, 1390 


A COMPARATIVE STUDY OF STRESS - FIELDS OF BLADE - DISC 
AHACHMENT CONFIGURATIONS OF A TURBINE BLADE 


A Thesis Submitted 
in Partial Fulfilment of the Requirements 
for the Degree of 
MASTER OF TECHNOLOGY 


by 

PITAMBAR CHAUDHARY 


to the 

DEPARTMENT OF MECHANICAL ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY, KANPUR 

OCTOBER, 1990 



Asti 

J_ ■ ^ 

I^i), 't I 



C 





CERTIFICATE ■ ' 

This is to certify that the thesis entitled "A Conparatlve Stw 
of Stress-Fields of Blade-Disc Attachment Configurations of 
Turbine Blade by Pitambar Chaudhary is a record of work carri 
under our supervision and has not been submitted elsewhere for 
degree. 


Dr N.S.Vyas" ''' 
(Assistant Professor) 
Dept, of Mechanical Engg. 
I . I . T . Kanpur 


Dr. P.M. Dixit 
(Assistant Professor) 
Dept, of Mechanical Engg. 
I I.T. ’Kanpur 


October, 1990 



ACKNOWLEDOEhEINT 


I wish to express my deep sense of gratitude and 
indebtness to Dr. P.M. Dixit and Dr. N.S. Vyas for their 

inspiring guidance, invaluable suggestions and constructive 

criticism. They were always a constant source of 
encouragement throughout my thesis work. 

I am thankful to Dr. N.N. Kish ore for allowing me to use 

the P.C. facilities in ESA Lab. 

I thank Mr. A.V. Rao, Mr. A.K. Keshari and Mr. K. Suresh Babu 
for their help rendered in completing my thesis. I feel 
happy to recall my association with all my friends 

during my stay at I.I.T. Kanpur. 

— Pitambar Chaudhary. 



CONTENTS 

Page 

LIST OF FIGURES ii 

LIST OF SYMBOLS iv 

ABSTRACT vii 

CHAPTER I INTRODUCTION 1 

1.1 Root region of a turbine blade 1 

1.1.1 Geometry of the root region 1 

1.1.2 Stresses in the root region 3 

1.2 Literature survey on the contact problem 4 

1.2.1 References using finite element 

method 4 

1.3 Objective and scope of the present work 6 

1.4 Plan of the thesis 13 

CHAPTER II THE FORMULATION - 14 

2.1 The blade 14 

2.2 The disc 14 

2.3 FEM Formulation of the contact problem 17 

2.3.1 FEM Solution procedure 17 

2.3.2 FEM Formulation of the plane stress 18 

2.3.3 Formulation of the contact problem 23 

2.3.3. 1 Contact condition 25 

2. 3. 3. 2 Application of contact 

condition 27 

2. 3. 3. 3 Substructuring 30 

2. 3. 3. 4 contact conditions 

application procedure 32 

2.3.4 Evaluation of stresses 33 

2.3.4. 1 Stresses in the contact 

region 33 

2. 3. 4. 2 Stresses at the corner point 34 

CHAPTER III Results and discussion ^ 37 

3.1 Stresses in the blade and disc 38 

3.2 Analysis of the root region 38 

3.2.1 Results for the tee root 41 

3.2.2 Results for the straddle root 44 

3.2.3 Comparision of the tee and strdlle 

root 49 

3.2.4 Contact stresses and displacements 56 

CHAPTER IV Conclusions and 'suggestions for the future work 70 

f . 


References 


71 



(U) 


LIST OF FIOURFS 


FIGURE 

NUMBER 

TITLE 

PAGE 

1 . 1 

The Tee-root joint between the blade and 



the disc 

7 

1.2 

The straddle root joint between the blade 



and the disc 

8 

1.3 

Root domain for the contact problem 



(Tee-root > 

11 

2. 1 

Blade mounted on a rotating disc 

15 

2.2 

Body in plane stress 

19 

2.3 

Discretisation of the domain 

19 

2.4 

Two bodies A and B in contact 

24 

3. 1 

The Tee-root joint between the blade and 



the disc 

39 

3.2 

Variation of the radial and circumferential 



stresses in the disc with radial co-ordinate 

40 

3.3 

Finite element mesh for tee-root joint 

42 

3.4 

Variation of the equivalent stress at the 
corner with respect to thickness to length 



ratio 

43 

3.5 

Straddle root joint between the blade and 



the disc 

46 

3.6 

Finite element mesh for t.he straddle root 



joint 

47-48 

3.7-3.11 

Variation of the equivalent stresses at the 

two corners with respect to the ratio of 



the contact shoulder lengths 

50-54 

3.12 

Variation of the equivalent stress in the 
optimum straddle root with respect to T/L 



ratio 

55 

3. 13 

Variation of non-dimensional equivalent 

stress in the optimum straddle root with 




( iii ) 

respect to T/L 57 

3.14 Variation of the normal displacement along 

the contact length (for tee-root) 59 

3.15 Variation of the tangential displacement 

along the contact length (for tee-root) 59 

3.16 Variation of the normal stress along the 

contact length (for tee-root) 60 

3.17 Variation of the shear stress along the 

contact length (for tee-root) 60 

3.18 Variation of the normal displacement along the 

contact length (for straddle root, model I) 61 

3.19 Variation of the tangential displacement along 

the contact length (for straddle root, 
model 1) 62 

3.20 Variation of the normal stress along the 

contact length (for straddle root, model I) 63 

3.21 Variation of the shear stress along the 

contact length (for straddle root, model I) 64 

3.22 Variation of the normal displacement along the 

contact length (for straddle root, model II) 65 

3.23 Variation of the tangential displacement 

along the contact length ( for straddle 

root, model II) 66 

3.24 Variation of the normal stress along the 

contact length (for straddle root, model II) 67 

3.25 Variation of the shear stress along the 

contact length ( for ^straddle root, 

model ID 68 



(IV ) 


LIST OF SYMBOLS 


A 

Area 

of 

cross 

-section 

A' ,B' 

Corner 

po 

ints 


B 

Width 

of 

the 

neck in 

bi^bg 

Width 

of 

the 

necks in 

C 

Thi ckn 

ess 

of 

the disc 


Thi ckn 

ess 

of 

the disc 


joint 





Constants 



E 

Young ' 

s modulus of ela 

F 

2 

Axial 

force 


h 

Axial 

thi ckness 

L 

Length 

of 

the 

shoulder 

lirlg 

Length 

of 

the 

shoulder 

lb 

Length 

of 

the 

blade 


x,y,z 

Z 


tee-root 
straddle root 
arm in tee-root 
arms in the straddle root 


in 


tee-root 

straddle 


root 



Outer 

radius 

of the disc 



Inner 

radius 

of the disc 


r,d 

Polar 

co-ordinates 


T 

Thi ckn 

ess of 

the shoulder 

in 

t^,tg 

Thickness of 

the shoulders 

in 

u, w 

Displacements 

X and Z dire 

cti 


tee-root 

straddle 


root 


Cartesian co-ordinates 

Distance along the blade length 


a 

Angle made by 

tangential 

di rect ion 

at 


contact node 

with X-axis 



r & 

Radial and circumferential 

strain 


P 

Density of the blade and 

disc material 


Angular veloc 

ity 




Co-efficient 

of friction 




Clearance 




X 2 

Clearance in 

X and Z-direclions in 

tee-root 

x1' x2 

Clearance in 

X-direction in 

straddle 

root 

21" z2 

Clearance in 

Z-direction in 

straddle 

root 

r 

System boundary 




the 



(V ) 


Boundary on which forces are specified 

O' , 0 ^ Radial and circumferential stress 

r & 

O' ,T Normal and shear stress 

n 

O' Stress at the Gauss-integration point 

g 

o Stress at the nodal point 

c 

O' Equivalent stress at the corner point 

O' Non-dimensional equivalent stress in the optimum 

eq 

straddle root 

C(o') Equivalent stress in the optimum straddle root 

opt s ^ 

(O' )_ Equivalent stress in the tee-root 

eq T ^ 

CB3 Matrix containing the derivative of shape 

functions 

CD3 Matrix containing material constants 

CK3 Stiffness matrix 

CN3 Matrix containing shape functions 

fB^> Body force vector 

Cf> Traction force vector 

fF> Force vector 

{U> Displacement vector 

io} Strain vector 

■Co'> Stress vector 



SUBSCRIPT 

A 

B 

e 

i » J » 
n , s 
T 


Cvi) 


AND SUPERSCRIPT 

Body A 
Body B 
Elemental 
Dummy indices 

Normal and tangential direction 
T ranspose 



ABSTRACT 


This work is a comparative study of stress-fields of the 
two most common blade-disc attachments of a turbine blade, 
namely tee and straddle. The study is restricted to only 
steady case i.e. it takes care of only the rotation of 
the disc but not the vibrations of the blade. For 
simplicity the blade is considered parallel to the disc 
and the placing of the root is assumed to be axial. The 
blade-disc attachment region is analysed as a plane stress 
contact problem by finite element method. The individual 
analysis of the blade and the disc provide the boundary 
conditions for the contact problem. 

A finite element program is developed to find the 
equivalent stress at the corner points as well as stresses 
and displacements in the contact region. First the 


straddle 

root is optimised 

so 

that 

the equivalent 

stress 

at 

the 

corner points of 

both 

the 

shoulders become 

equal . 

Then 

the 

optimum straddle 

root 

is 

compared with 

tee-root 

of 

the 

same thickness to 

length (T/L) ratio for 

various 


values of this ratio. Finally normal and tangential 
stresses and displacements are presented for typical cases 
of tee and straddle roots. 
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CHAPTER I 

INTRODUCTION 

11 ROOT REGION OF A TURBINE KJVDE ; 

The blades in a turbomachine are connected to the 
rotor by an attachment, called as "root region". The 
reliability of a turbomachine is affected by the 
reliability of its blades as well as the root region, as 

the failure of any of the blades or the root region 
leads to the failure of the turbine. Hence, a turbine 
blade and its root region must be strong enough to 

withstand different type of forces acting on it. 

11.1 GEOMETRY OF THE ROOT REGION : 

There are various ways of fixing turbine blades 

to the rotor. Any of the following methods can be used 

depending upon the application: 

(i> Rotor and blades are machined as a single piece. 

Since this is an expensive arrangement, it is used in 
very special cases where the speed and working temperature 

of the rotor are very high. 

s. 

(2) Blades are welded to the rotor. This arrangement is 

used in some gas turbines. Like the first type this 

arrangement also has light weight. But it is 

unsatisfactory from the maintenance point of view. 

(3) Blades are inserted radially in a circumferential 

groove in the rotor. This arrangement is normally used in 

impulse turbine or in a low pressure stage of large 

turbines. In this case common shape of the roots are 
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dove-tai 1 , tee 

or straddle. 




(4) Blades are 

inserted 

axial ly 

in the axial 

slots 

spaced 

uniformly around 

the 

periphery 

of the 

rotor . 

This 

arrangement is 

normally 

used for 

fixing long 

blades 

in a 

low pressure stage of 

multistage 

turbine . 

The 

common 

shapes of the 

root for 

this case are fir-tree and. 

ball 

and shank type. 







For the radial placing of the blade roots the 


circumferential groove is widened at one place. This is 


called 

the gate. 

The blades are inserted one 

at a 

time 

through 

the gate 

and slid around the slot 

until 

the 

circumference is 

nearly full. Finally the gate 

is 

closed 

by some 

special 

stop-piece. For balancing reason 

another 

similar 

gate is 

made at 180° apart. Since 

the 

contact 

surface 

for the 

dove-tail root is inclined. 

the 

joint 


becomes loose when the turbine is subjected to repeated 
heating and cooling. On the other hand shoulder of tee 
and straddle roots give better grip to the blade. The 

straddle root is preferred when additional contact is 
required. 


For axial placing of the blade roots the 
periphery is divided into as many points as the number 
of the blades in a single row. Sufficient axial spacing 
between the adjacent rotors is necessary to facilitate 
the insertion of the blades. In the case of large 
centrifugal force the axial dimension of the root should 
be wide enough to provide sufficient contact area. 
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112 STREdSCd W THE ROOT REGION 

In the turbomachines the support provided by the 


disc 

to the 

blades is 

more like a 

rigid 

support. As 

a 

result 

maximum 

stresses 

occur in the 

root 

region of 

the 

blade . 

Further 

, because 

of the existance 

of the 

sharp 


corners, there is a stress concentration. Therefore, it is 
extremely important to carry out the stress analysis of 
the root region in general and critical corner points in 
the particular. 

The stresses induced in the turbine blade are 
generally of two kinds. Those arising due to the rotation 
of the disc do not vary with time. These stresses have 

been termed as “steady stresses*' in this thesis. The 
blades are also subjected to nozzle forces. These forces 

cause vibration of the blades and induce time dependent 
stresses in them. These can be called as "non-steady 

stresses". The complete dynamic analysis involves 

determination of both the steady as well as non-steady 

parts. The problem of the stress analysis of the root 
region, because of its complex geometry, is quite 
difficult even for the steady case. Addition of non-steady 
effect will make it even more difficult. Therfore, to 

keep the problem to the manageable level only steady-state 
stress analysis is considered. 

While analysing the root region of a turbine 
blade its interaction with the surrounding slot of the 

disc has to be taken into account. This part of the 
disc does not behave any more rigidly than the root of 
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the 

blade . 

Thus 

the problem 

of 

the root 

region 

has 

to 

be 

treated 

as a 

problem of 

the 

stress 

analysis 

of 

two 

elastic bod 

ies in 

contact. ! 

Such 

problems 

are 

called 

as 


"contact problems". 

12 LITERATURE SURVEY ON THE CONTACT PROBLEH : 

The first analytical approach to find the stresses 
and displacements in elastic bodies in contact with each 

other was proposed by Hertz HID. This was a three 
dimensional problem based on certain assumptions. The two 
most important assumptions were that the area of contact 
was very small in comparison to the dimension of the 
bodies and the contact was frictionless. The three 
dimensional Hertz contact problem was simplified to an 
integral equation by 1. Ya. Shtaerman C23. Using this 

approach several contact problems were solved by 
Mushketishvili C33 and Gladwell C43. However, all the 


contact 

problems 

so] 

Ived so 

far 

have simple 

geometries 

. No 

wonder , 

the probl 

erri 

of the 

root 

region of 

a turbine 

blade 

does nc 

»t figure 

in 

this 

list. 





12.1 REFERENCES USING FWrTE ELEMENT METHOD : 

There have been many attempts to use numerical 

methods to solve contact problems. The most widely used 
method is the finite element method <FEM). One of the 

earliest important reference as on use of FEW in solution 
of the contact problem is that of Francavilla and 
Zienkiewicz CS3. They solved the frictionless contact 
problem by treating it as quasi-linear problem. N. Okamoto 
and M. Nakazawa C6D proposed a finite element method based 
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on load incremental theory taking -friction into account. 
Hung and Sauxe C71 used FEH and mathematical programming 
techniques to solve frictionless contact problems. The 
problem was solved as a non-iterative optimisation problem 
where the contact conditions were linearised in a suitable 
fashion. Sachdeva and Ramakrishnan extended the method of 

Francavilla and Zienkiewicz to take care of force boundary 
conditions C8D and frictional effects under proportionate 
loading C9I1. Torstenfelt CIOD proposed a general purpose 

finite element computer program to solve contact problems 
with friction. Rahman et. al . C113 used finite element 

incremental procedure for contact analysis where one of 
the bodies is rigid. The contact condition at a node is 

established when it penetrates or touches the rigid body. 
Bathe and Chaudhary C12D presented a procedure for the 

analysis of plannar and axisymmetric contact problems 

involving sticking, frictional sliding and separation under 
large deformation. Contact of a cylinder with a rigid 

plane, motion of a rubber sheet in a converging channel 
and analysis of a buried pipe were solved as 
illustrative problems. Chandrasekhra et. al. C13D solved the 

contact problem by imposing geometric constraints on the 
pseudo equilibrium configuration defined as a configuration 
at which the compatibility conditions are violated. 


The first reference on use of FEM in analysis of 


a turbine blade 

root 

region 

is 

that 

of 

Cham and 

Tuba C143. 

They 

have used 

triangular 

elements 

to 

study the 

fir-tree 

root . 

Their results 

show 

that 

the 

coefficient of 

friction 



has only a marginal effect on stresses and displacements 
in the root region but the effect of the clearance 
between blade root shoulder and disc arm is quite 
significant. However, the detailed results are not 
presented. Srivastava C15D analysed the tee— root to study 
the effects of rotor speed and the clearance around the 


root 

shoulder on 

the contact stresses 

and 

displacements . 

The 

four noded 

isoparametric 

elements 

were 

used. 

13 OBJECTIVE AND 

SCOPE OF THE 

PRESENT 

WORK : 



The present work is an attempt to carry out the 

steady-state stress analysis of the blade-disc attachment 
configuration of a turbine blade. The study is restricted 
to only two configurations namely tee-shaped and straddle 
roots. The objective is to compare these configurations 
from the point of equivalent stress at the corner points. 
The configurations are shown in Figs. 1.1' and . The 


straddle 

has 

two necks. 

We 


have considered 

two 

different 

models of straddle 

root . 


In the first model 

the 

width of 

both 

the necks are 

equal 

while in 

the second 

the width 

of 

the lower neck 

is 

s, 

reduced 

appropriately. 

This has 

been 

done keeping 

in mind 

the fact 

that 

the 

lower neck 

is 

subjected to 

less 

axial load 

compared 

to 

the top neck. 

The straddle 

root 

also has two shoulders 

and hence 

two 

critical points. The 

equivalent 

stresses 

at 

these two 

corner points depend 

on 

relative 

lengths 

and 

thicknesses 

of 

the shoulders. 

In 

the optimum 

design 

of 

a straddle 

root, equivalent 

stress 

at both 

the corner 


points should be equal so as to make effective use of 




For the blode 

B = Width of the n^ck 
L = Length of the shoulder 
T = Thickness of the should 
Clearance in x direction 
Clearance in z direction 

For disc 

C = Thickness of disc arm 


AB ond CDEF boundaries define the root domain. 
Inner radius of disc = 50 mm. 


The TeTfoot joint between the blade and disc 





AB ond CDEF boundories define the blode ond disc straddle 


b ■= Width of the 
I = length of the sholder 
t = Thickness of the should 
€x= Clearance in x direction 
€z= Clearance in z direction 


Subscript 1 to the upper 
shoulder 

Subscript 2 to the lower 
shoulder 

F 



Fig 1.2 Straddle root joint between the blade 
and disc 



3 


shou 

ilders . 

The 

optimum case 

IS obtained 

by 

the 

ratio 

of 

lengths of the 

shoulders . 

The 

root 

with 

the 

optimum ratio of 

shoulders 

length 


is henceforth referred as "optimum straddle". This optimum 
straddle root is compared with a tee-root. Since there 

has to be some common ground for comparison, the total 

contact area in both the roots is kept the same. This 
means the sum of the lengths of the two shoulders of 
the straddle root is made equal to the length of the 

tee-root. The thicknesses of the two shoulders of the 

straddle root are taken in proportion to their lengths 
and their sum is made equal to the thickness of the 

shoulder of the tee-root. Similar policy is adopted about 
various clearances. Thus the relations between 
geometric parameters of straddle and tee root are 
by 


various 

given 


+ Ig = L 
t, + tg = T 


where 


« + : 


where 


^ + lit SS 4P 

Zl ^2 == 


( 1 . 1 ) 


where 
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( see 

the Figs. 1.1(a) and 

1.1(b) ). 

The 

comparison 

of 

the 

tee and straddle roots 

is made 

for 

various values 

of 

the 

thickness to length (T/L) ratio. 





The blade is normally staggered with respect to 


the disc. This makes the actual problem three-dimensional. 
However to keep things simple, we assume the blade to be 
in the plane of the disc. Further it is assumed that 
the blade is inserted axially. Because of these two 
assumptions, the problem can now be considered as a plane 
stress problem. Finally the blade is idealised as a bar 
of uniform rectangular cross-section and the disc thickness 
is assumed to be uniform. First the blade is analysed as 
a bar subjected to variable body force resulting from the 
rotation of the disc and the disc is analysed as a body 
in the state of axisymmetric plane state and subjected to 
variable body force in the radial direction. But this 
analysis becomes invalid in the root region because of the 


change in 

geometry 

and 

interaction 

between 

the blade 

and 

disc. To 

analyse the root region. 

a domain 

consisting 

of 

appropriate 

portions 

of 

the blade 

and disc 

surrounding 

the 

actual contact area 

is 

chosen as 

^shown in 

Fig. 1.2 . 

The 

boundaries 

of the 

domain 

adjoining 

the blade and the 

disc 

(i.e. the 

boundaries 

AB, 

CD, DE and 1 

EF) have 

to be at 

a 


sufficiently large distance from the contact area so that 
the analysis of the blade and the disc, as done above, 
remain valid upto these boundaries. In fact, the stress 
boundary conditions on these boundaries are provided by 
the blade and the disc analysis. The root region is 



z 



Fig.j.3 Root domain for the contact problem. 

(Tee root) 
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analysed as a plane stress contact problem using Finite 
Element Method. 


A computer program is developed to determine the 


equivalent stress at 

the corner po 

ints 

of the 

shoul 

der as 

well 

as to evaluate 

the stresses 

and 

displacements 

in the 

conta 

cl region. First 

the tee-root 

i s 

analysed 

to 

obtain 

the 

variation of the 

equivalent st 

ress 

at 

the 

corner 

point 

with respect to 

thickness to 

length (T/L) 

ratio 

. Then 

for 

a given T/L rati 

0 , the models 

of 

straddle 

root 

are 


analysed to find the optimum ratio of shoulder lengths 
(1^/lg). At the optimum value, the equivalent stresses at 
both the corners becomes equal. Finally the equivalent 
stress in the optimum straddle root is compared with the 


equivalent stress in 

t h© 

tee-root 

with 

the 

same 

T/L 

ratio. This is done 

to 

obtain the 

range 

of 

T/L 

ratio 


for which the straddle root is better than the tee-root 
i.e. has a lesser value of the equivalent stress. Variation 

of the contact stresses and displacements for both the tee 
and straddle roots is reported for the lowest T/L ratio 

from the favourable range. 

s 

This study has several limitations in the form of 

approKimations and idealisations made as under. 

(1) The study is restricted to only steady conditions. 

Under unsteady conditions the range of T/L ratios for 

which the straddle root becomes better than the tee-root 
may change. 

(2) The blades with tee or straddle roots are usually 
fixed radially. But this type of attachment makes the 
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problem three-dimensional. To reduce the problem to a 
state of plane stress the blades are assumed to be fixed 
axially. 

(3) The blade is analysed as a rectangular bar of 

uniform cross-section under axial loading. But the actual 

blade is tapered, twisted, and has an aerofoil 
cross-section. 

(4) The disc is considered as an axisymmetric body of 

uniform thickness. However, in practice, the disc with 

axial slots is not axisymmetric. Moreover, the thickness 

of the disc varies in the radial direction. 

14 - PLAN OF THE THEdlS : 

The formulation of the blade, disc and contact 
problem are discussed in Chapter II. The solution procedure 

for the contact problem is also described here. The 

Chapter III contains the results and discussions. The 

conclusions and suggestions for future work are given in 


Chapter IV. 
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CHAPTER II 

THE FORMULATION 


2.1 THE BLADE : 

The blade is treated as a bar of uniform 
rectangular cross-section mounted on a f'otating disc, 
(Fig. 2.1). The axial stress induced in the blade due to 
disc rotation is evaluated as follow. 

The centrifugal force on a small element of length 
dZ at a distance Z from the root of the blade (Fig. 2.1) 
IS given by 


dF = p 


A (R,+ Z) dZ 
1 


Z IS given by 

1 . 




force 

on 

a 

(R.,+ 

Z) dZ 


R 7 4* 




2 

iz=z 

a 

on 

this 


distance 


— f> 




(1. - Z ) + 
b 


(iVz^) 

b 


( 2 . 1 ) 


2£ THE DISC : 

The disc is assumed to be of uniform thickness 
clamped at its inner boundary on a shaft which is rotating 
at a constant angular velocity about the Y-axis (Fig. 2.1). 
The disc is in a state of axisymmetric as well as plane 
stress. There is only one nontrivial equation of motion. In 
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Fig. 2.1 Blade mounted as a rotating disc. 


\6 


terms of stress components, it is given by 


d & 


''a. r & ^ z 

— + + f> r — 0 


(E.E) 


dr r 

Since it is one equation in two unknowns, it 

cannot be solved in the present state. However, it can be 

solved by expressing it in terms of the radial displacements. 

For the axisymmetric plane stress problem, the 

stress-strain relations are 

E 


1 - V 




= 


1 - 


( + V £ ) 

e r 


(2.3) 


and the strain displacement relations are 

d u u 

r r 

^ = • 4} = 

r ' 0 


(2.4) 


dr r 

Substituting the relations (2.3) and (2.4) into (2.2) we get 


1-fc* 




(ru)| + ^< 1 ^ r = 0 


the solution of this equation is 

1 - 1 *^ 2 ^ 

u = - p ^ C. r + — - 

r r- 8 1 r 


(2.5) 


( 2 . 6 ) 


where and Cg are constants of integration. 

Substituting the expression (2.6) for the radial 

displacement into the equation (2.4) and (2.3) we get the 


following expression for stresses 


<» = C. 

r 1 


1-s* 


1-Hi 


3+fc» 


8 


2 2 
p uy r 


(2.7) 


"a = <^1 


1-0* 


'2 


1+V> 


1+31* 


8 


2 2 
p iJy r 


( 2 . 8 ) 
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The constants and Cg are to be deteriraned from 

the boundary conditions. Since the shaft to which the disc 

is attached is relatively more rigid than the disc, at the 

boundary r = Rg the boundary condition becomes 

u^ = 0 at r = Rg . (E.9> 

The boundary r = R^ is stress free. Therefore 

= 0 at r = R^ (£.10) 

Using these conditions, we get 

C = 1- + p R (2.11) 

1+v R^ BE 

1 

with 

R^ R^ - (3-K*)R^ - (1+v)R^ , 

C = - 1 ^ J ^ (2.12) 

8E *■ (1+w)R^ + (1-w)R^ J 

2.3 FEM FORI-RAJiiTION OF THE CONTACT PROBLEM : 

As stated earlier the root region of the turbine 


blade 

is analysed as 

a contact 

problem. 

The domain 

of 

the 

problem is shown 

in Fig. 1 

.2. The 

boundary conditions 

for 

the problem come 

from the 

analysis 

of section 

2.1 


and 2.2. Since the domain is a plate like region and 

stresses on the boundaries are constant across the 
thickness, this contact problem can be treated as a plane 
stress problem. 

2 . 3.1 FEM SOLUTION PROCEDURE : 

In FEM formulation, the governing equations (partial 
differential equations) are converted into an equivalent set 
of algebraic equations. There are two common procedures 
available for converting the governing equations into 
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equivalent algebraic equations. These are (i) the variational 
procedure and (ii) the Galerkin procedure (or the method of 


weighted residues). The 

solution domain 

IS 

divided 

into 

several 

small elements 

of chosen shape 

and 

the 

unknown 

field 

variables are 

expressed in 

terms 

of 

assumed 

approximating functions 

(or interpolation functions) 

within 

each 

element. The 

approximating functions 

are 

defined 

within 

each element in 

terms of the 

field 

variables at 

specified points called 

as nodes The 

nodal 

values of the 

field 

variables and the interpolation 

functions for 

the 

element 

completely define the behavior of 

the 

field 


variables within the element. For the finite element 

formulation of the problem, the nodal values of the field 
variables become the unknowns. The algebraic equations of 
the element involve integration of the interpolation 
functions or their derivatives over the element. All such 

elemental integrals are assembled to obtain the final matriK 

equation. The unknown nodal variables are calculated by 

solving the assembled matrix equation for the whole 
solution domain. 

2.3J2 FEM FORMULATION FOR PLANE STRESS ; 



Consider 

a body 

O of 

uniform 

thickness h 

bounded 

by two 

parallel 

planes 

Y = -h/2 

and 

Y 

= h/2 and a 

closed 

boundary 

r (Fig. 

2 » 2 ) • The part 

of 

the 

boundary on 

whi ch 

external 

stresses 

act is 

denoted 

by 

q • 

The body 

is in a 

state of 

plane 

stress . 

Although 

y 

is 

not zero. 

it is 

not required 

in the analysis 

• 

Only the 

strain 



1 3 



Fig. 2.2 Body in plane stress 



Fig. 2. 3 Discretisation of the domain. 
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components 


in 


X-Z 


plane are 


needed . 


The 


strain-displacements 

du 

* - 1 — r 


relation 


in this plane are 


+ 5^ 


xz 


In matrix form 

<«>=CL3{u> 

Where £«> = ■[« « y 

X z xz 

fu> = fu w)"*^ 


C L 3 



( 2 . 13 ) 


( 2 . 14 ) 


Then stress strain relationship for a homogeneous, 
isotropic, linearly elastic material under plane stress 
conditions is given by 

C«y>=CD3<:«> (2.15) 

Where C > = C 


and 


C D 3 = 


1-v 

The strain energy is 


^ V 0 

V 1 
0 0 - 


V 1 0 

v »— 1 
2 


(2.16) 


U=:J|{*>''^{4»JhdA 


(2.17) 


If 




and 


{••i-C:! 


(2.18) 


represent the applied traction on r and the body force 



2.1 


(per unit volume) respectively, then the work done by 
the external forces is given as 

Wp=r fu>^<B^>hdA + r {u>^(:f>hdl <E.19) 

■'r- 

q 

The expression for Potential Energy of the system is given 
as 

= U - (2.E0) 

p p 

Using Eqs. (2.15), (2.17), and (2.19) 



C ic }^C D J{ ^ >hdA 


u >hdA 

■*A 


Jf u f >hdl 

r 

q 

( 2 . 21 ) 


Let the domain be divided into N number of 

8 

elements (Fig. 2. 3). Over a typical element e, let 

r e 1 

*^1 

N. 0 . . . . N 01 w® 

1 nil 

0 .... 0 N J ' -e ^ 

e 
w 

n 

I J 

= ^ “ >no<lal 

be the approximation for the displacement vector, where n is 

the number of nodes per element u^®, w^®, u^^®, w^® 

are the displacements at the nodes and are the 

shape functions in the form of polynomials of the 

coordinates x and z. Then from equation (2.13), we get 

f*>=CB3<u> (2.23) 

nodal 

Where the matrix C B 3 contains the derivatives of the 

shape functions. Substituting Equations (2.22) and (2.23) in 


C u > = 



to the expression (S.S1) for the pot.ential energy, we get 
N 


n =) 4'‘'f 

e=1 e 


podal ® ^ C D a c B 3 { u dfl 


N 


N. 


LK " 


. ,C N D'f B. >dA 
nodal f 


^ A 

e=1 e 


” ^ ^nodal 

e=1 * o 


C N 3 C f > 


dl 


(2.E4) 


Here is the boundary of an element which is coirimon 

with P and N. is the number of such elements 
q D 

(Fig. 2.3). 

Setting the first variation of to zero, we get 


CK3fU> = -CF> 


(2.25) 


Where 


N 


CK3=^CK3' 


e=1 


(2.26) 


K 3®= J C B 3^ 


CD3CB3hdA 


(2.27) 


{ F > = 


V C F 
e=1 


(2.28) 


{ F >®= Jc N 3^C >h dA + J C N 3 ^c f >h dl 


(2.29) 


A 


and f U > is the vector of nodal displacements for the 


entire domain 
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In the equation <2.26) and (2.28), it is understood 
that the sum is performed after ejtpanding C K and { f }® 

to the full size. 

2 . 3.3 FORMULATION OF ThC CONTACT PROBLEM : 

Two bodies A and B which are in the state of 
plane stress are in contact with each other as shown in 
Fig 2.4. The finite element equations for these bodies are 

[•<«]{ }'{ 

[ S ] { ‘'b } = { >^6 } ‘2.30b) 

where [ ] , [ Kg ] • { } . { Ug } ' { 

stiffness matrices, displacement vectors and force vectors 

of bodies A and B respectively. The force vector includes 
the contact forces also. 

The equations (2.30) cannot be solved as at the 
contact nodes, neither displacements nor forces are known. 

To overcome this difficulty, contact conditions are used. 

Before applying the contact conditions, it is 
convenient to combine the equations (2.30a) and (2.30b) 
into one single equation. 



or 


C K U > 


{ F > 


(2.31) 





£. 3 . 3.1 CONTACT CONDITIONS : 

There are two sets o 
the unl’nown contact forces 
displacements. Let node of 

the node of body B (Fig. 

respectively the normal and the 
contact node. Then the condit 
Newton’s third law, are 


25 


contact conditions. 

one 

on 

and 

the 

other 

on 

the 

body 

A be 

in conta 

ct 

with 

.4 ). 

Let n 

and s 

denote 

tangen 

tial 

di rections 

at 

the 


ions on forces, as given by 


( ""a ] . ' - ( '"b ) . <2-32) 

'■ ■'in '■ jn 

( '"a ] . “ ■ [ '"b ) . <2.33) 

^ -'jS 

The compatibility of displacement components gives 
the second set of condition 


^ in ^ jn 

where « denotes the clearance 

The compatibility condition 
direction is valid only when there 
contact nodes. The slip depends on 
normal force at the node and the 


(2.34) 


(2.35) 


■for the tangential 
is no slip between the 
the magnitude of the 
coefficient of friction 

slip, the 
should be 


between the bodies in contact. In the case of 
condition (Eq. 2.35) on tangential displacement 
replaced by the following relation 



2 € 

I ( '"bJjJ I [ '2-36) 

This condition can be applied to the body A also. 

The contact conditions <2.3S)-(2.36) are in terms of 
the components with respect to normal and tangential 

directions, while the finite elements equation (2.31) are 
with respect to X-Z coordinate system. Therefore, first, 

these conditions must be expressed in terms of the X-Z 


components. From 

the 

Fig. 2. 4, we 

get the 

following 

relations between 

the 

n-s components 

and X-Z 

components . 

(F.). = 

A in 


. ^ Sin « + (F. ) . 
IX A 12 

Cos « 

(2.37) 

(F-). = 

A IS 

^^A> 

. Cos a + (F- ) . 

IX A 12 

Sin a 

(2.38) 

Here a is 

the 

angle which s- 

■directions 

makes with 


x-axis. Similarly, the normal and tangential components of 

other forces and displacements can be expressed in terms of 

the X-Z components. Substituting the Eq. (2.37)-(2.38) and 

similar relation for = » . 

B jn B js A in A is B jn 

and into the Eq. (2.32)— (2.36) , we get the following 

conditions 

-(F^)i^ Sin « + <F^).^ Cos « - <Fg) Sin « -h Cos « = 0 

(2.39) 

(F^)i^ Cos « + (F^).^ Sin « + (Fg) Cos « + (Fg) Sin « = 0 

(2.40) 

(U. ) . Sin « - (U_ ) . Cos « - (U„) . Sin « +(U_) . Cos « = jp 
A IX A 12 B jx B j 2 

(2.41) 

(U.). Cos « + <U-). Sin <* - (Ujj). Cos a - (U^) . Sin « = 0 

H IX H IZ » JX D JZ 


(2.42) 
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(ss cosot +p sn Sinot )(F ) + (ss sinoi - fj sn cosotXFr,). = 0 

B jx B j: 

<S.43) 


where 

ss is sign 

of 


and 

sn is sign of 

(F 

’b 

that 

in case of 

no 

slip 

the 

condition 

(2.42) 

applies while 

the condition (2. 

43) 

is to 

be 

used when 

there 

is 

slip. 

2-3.3£ 

APPLICATION 

OF 

CONTACT 

CONDITIONS ; 





II can be seen that cor responding to every contact 
pair (irj), there are four equations in the set (S.31) : two 
equations - (Hi-1)th Eq. for x-direction and <2i>th Eq. for 
z-direction, for the body A and the remaining two (2j-1)th 
Eq. and <2j>th, for the body B. At stated earlier, neither 
the displacement associated with these nodes nor the forces 
(i.e. the right hand side of these equations) are known. 


Therefore , 

these equations 

can 

not be 

solved 

as such. 

However for 

every contact 

pair 

, there 

are 

four contact 

conditions. 

Our strategy is 

to 

replace 

these 

equations with 


unknown right hand sides by the contact conditions which 

have the known right hand side. The matrix. operations 
required for this process are described below. 

First we replace the (2i-1>th Eq. of body A by 

the Eq.<2.39). Since and a*"® equal to the 

left hand sides of <2i-1)th and (2i)th Eq. of the body A 
and (Fjj) • and <Fo) 4 , are equal to the left hand side of 

<2j-1)th and (2j)th Eq. of body B, the condition states 


that 



2 ^ 


(-sin«) (2i-1)th Eq. of body A 
+ (cosot) (2i)th Eq. of body A 
+ (-sirKx) (2j-1>th Eq. of body B 
+(cosa) (2j)lh Eq. of body B =0 
Thus after replacement, the (Ei~1)th row of the CK.D becomes 
equal to 

<-sino« ) (2i-1 > th row of 


+ ( cosoi ) (2i ) 

th 

row 

of 

CKa^ 

+ (-sinoi ) (2j 

-1 ) 

th 

row 

of 

+ ( cosot ) (2j ) 

th 

row 

of 

CKgD 

The (2i-1)th 

row 

of 

CFa3, 

of 


^B- 


Next, the <2ilth Eq . of body A is replaced by the 
Eq. (2.40). The matrix operations for this replacement are 
similar. Now, the (2i) th row of is replaced by the 

linear combination of (2i-1)th and (2i)th rows of sind 

(2j— 1)th and (Ej) th rows of CKg3, the coefficients of a 
linear combination given by the Eq. (2.40). The (2i)th row 
of O'f course, becomes zero. Note that the Eq.(2.40) 

: tKia tr nttt Pk 

'A' 


involves the original (2i-1)th row of 


and not the 


modified one. Thus (2i-1)th row of ^isis to be stored 

before applying the first contact condition. 

The third condition is applied in the end. The 

reasons for this change of order will be explained later. 

So now, consider the fourth contact condition. Me use this 
to replace <Sj-1)th Eq. of body B. In case of no slip 
the condition is given by the Eq. (E.42). Since this Eq. is 
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in terms of displacement components, we don’t have to do 

any rove operation like the first two conditions while using 
it for 'eplacement. Instead, we simply have to replace the 
<2j-1)th row of zeros except in the columns 

(Ei-1), Hi), (Hj-1) and (2j). In these columns, the new 

coefficients will be respectively cosot, sinoi, ( — cosoi) and 

(-sina). The (Ej-1)th row of of course, becomes zero. 

In case if slip, the condition is given by the Eq. (2.43). This 
condition is similar to the first two conditions except 
that It involves only two Eq. : <2j-1)th and (2j) th 

equations of body B. So here, (2j-1)th row of 
replaced by the linear combinations of (Ej-1)th and (2j) th 

rows of CKg3, the coefficients of the linear combination 

given b> Eq. (2.43). The (2j-1)th row of CFgD is, of course, 

made ^itr >. Finally, the (Ej) th Eq. of body B is replaced 
by the ‘-bird contact condition (Eq. 2.41). Since it is in 
terms of the displacement components, we simply have to 

replace -he (2j)th row of ^y zeros, except in the 

columns (Ei-1), (Ei), (2j-1) and (2j) where the new 

coefficients will be sina, (-cosa), ( -sina) and cosa. The 

(2j)th row of ClFgT is made equal to * (the clearance). 

If we apply the third Eq. before the fourth, then 

the (2j)th Eq. of body B gets modified. In case of slip, 

the original {2j)th Eq. is needed to apply the fourth 
condition. So either we store the (2j)th row of EIKg3 and 

then apply the third and fourth condition in that order or 
we reverse the order and don't store the (2j)th row of 
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CHgD. The later alternative 
followed here. 


is more convenient 


and 


2.3.3 3 SUB-STRUCTURINQ : 

Since the number of contact nodes is initially 
unknown, an iterative procedure needs to be applied to the 
Eq.(2.31). Further imposition of the contact conditions 
renders the stiffness matrix unsymmetric. An iterative 

scheme involving large unsymmetric matrices r e quires a large 
storage and a considerable amount of computing time. To 
reduce both storage and computational time, the stiffness 
matrix is stored in skyline form and then decomposed to a 
much smaller size using substructuring. 

First step in substructuring is to rearrange the 
rows and columns of the global stiffness matrix as shown 
below 



Here the suffix I denotes the degree of freedom associated 
with contact nodes and the nodes at which non zero forces 
are specified, II denotes the degree of freedom associated 
with the free surface and internal nodes, while III denotes 
the degree of freedom associated with the nodes at which 



displacements are specified. If 


the prescribed displacements 


are 

zero 

then = 0- Further 

the 

nodal 

force 

at 


the 

free 

surf a 

ce and internal nodes 

are 

zero 

i . e . 

•^ii 

s 

0 . 

Thus 

the 

Eq. (2.44) reduces to 








'<1.1 

4 " "i.ii "'ll ' >^1 





(2 

.45) 



'‘11,1 **11,11 "11 ' 

0 




<2 

.46) 

El imi 

nat ing 

^II above 

equat 

ions we get 






^1,1 ^I ^ ’^I 





(2 

.47) 

where 

TTJ 

,I * ^ ^1,1 “ ^I,II ^II 

,II 

*^11,1 

) 


(2 

.48) 


To 

obtain the matrix 

*<1,1 

, K 

I, I 

is 

f 

irst 

extra 

cted 1 

n full form from the 

global sti 

f fness 

mat 

r ix 

in 


skyline form. Next K^.^. is obtained in skyline form from 

— 1 

the global stiffness matrix. Then the matrix K__ __ 

XXfxX XX7X 
t H 

IS obtained columnwise. If C C is its i column, then it 
is obtained by solving the system 

i i Column 


Thus 

the 

inverse 

of Kjj 

,II 

not explicitly 

found. Next 

the 

-1 

elements of the matrix K^. 

KlI,I 

are obtained 

by 

multiplying 

the rows of 

X y 

11 

the 

appropriate 

column 


*<11,1 • 

In this 

step as 

well 

as in the earlier 

step, we 

need the 

matrices 


and 

Kii j (these 

are 

transpose 

of each 

other) . 

But 

these 

matrices are 

not determined explicitly. 

Instead 

their 

rows 
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and columns are generated as and when required. Finally 


fn T IS obtained using the eq. (2.48). 

i IP -I 


Subsl r uc t ur ing 

thus leads 

to 

a 

smaller set of 

equations (2.47). This 

set is solved 

after 

applying the 

contact conditions. 

Di s placements 


are 

obtained by back 


substituting Uj in the quation 


U 


II 


*^ 11,11 ^ 11,1 


(2.50) 


Once all the displacements are known F,ii can 

be readily found using the third set of equation (2.44). 

23.3.4. CONTACT CONDITIONS AF1FH.ICATION PROCEDURE : 

The set of equations (2.47) are solved iteratively 

by applying the contact conditions in the following way 

(1) A set of nodal points is assumed to be in contact 

(it should include all possible contact points) and 
equations (2.47) are solved after applying the contact 
conditions where the fourth condition is the no slip 

condition, Eq.(2.42). 

(2) The normal nodal forces cannot be tensile. Hence all 
those nodes where normal nodal forces come out to be 
tensile are deleted from the possible contact zone in the 
next iteration. 

(3) The ratio of the tangential nodal force to the normal 

nodal force is calculated at all the node, in contact. If 
this ratio is greater than the coefficient of friction 

between the two surfaces at any node, it implies that 
this node is slipping in the tangential direction relative 
to the corresponding contact node. For the nodes whi 



slip. 


the contact condition (2.42) 
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IS replaced by (2.43). 


in the 

ne 

:;t 

1 terat 

ion . 

The 

other 

nodes are still 

gove 

rned 

by the 

cond 1 1 

ion (2 

.42) . 







Th 

e 

steps 

2 

and 

3 ar 

e repeated until 

all 

the 

norma 1 

for 

ces 

in 

the 

contact 

zone come out to 

be 

compres 

s 1 ve 

an 

d the 

rati 0 

of 

the 

tangential force 

to 

the 

nor ma 1 

for 

ce f 

or all 

the 

nodes 

in 

contact is 

ei 

ther 

equal 

or 

less 

than 

the 

coef f 

i c 1 ent 

of friction. 

Once 

the 

i te rat 1 

ons 

are 

over 

the 

stre 

sses 

are computed. 




2.3.4 EVALUATK^ Of STI^ISSES ; 


To 

evaluate the stress at a 

point , 

first 

the 

displacement 

vector C u } for the 

nodal 

element , 

to 

which 


the point belongs, is obtained from the solution vector 
{ U }. Then the stress vector is given by the equations 
(2.15) and (2.23). 


In 

sampling points as 
concerned. This is 
the inter element boundaries. 

2 . 3 . 4.1 STRESSES W THE CONTACT 
First, certain points 


are 

away 

from 

the 

nodes 

are 

Then 

the 

state 

of 

stress 

at 

by 

first 

considering 

1 them 

as 


stress vector € t > at that 


(2.51) 
the worse 
stress - is 
not continuous at 

REdON : 

on the contact surface which 
selected as sampling points, 
these points is determined 
a part of body A. The 
point is given by 


•C«'> = CD3IB>€u> 


nodal 


general the element nodes 
far as determination 
because stresses are 


are 

of 
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where <n> 


is the unit outward normal at that 


point to the boundary of body A. Similarly, the stress 
vector -C t } IS calculated by considering the point as the 
part of body B. Then the actual stress vector is 
determined by averaging it over bodies A and B. Finally, 
the normal and shear comp»onents of stress are calculated 
as follows 5 


€r 

n 


t n 

;< X 


t 


T 


n 


( 2 . 53 ) 



( 2 . 54 ) 


2.3.4-.2 STRESSES AT THE CORNER POWT : 

For the evaluation of stress at the corner points 
of the neck, (5auss integration points are selected as 

sampling points. The stresses at a node of an element 
are calculated by linear interpolation of the stresses at 
the Gauss points. The interpolation is done as follows. 
Similar to the approximation of the displacement, a stress 
component at a Gauss point (i^^) is approximated as a sum 
of the products of the shape functions (N^ ) with its 



nodal 


val uf" (O' ) 

c j 
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where i = 1 r 2 , -■ 


I = 1 , 2 , — 4 , 

= < 0 .). <N.) 


Thus at the four Gauss points 

(Ng)^ + (o- )^ (N„), 


c'3 '”31 


" <-cU 


• Jt 

By inverting this set of equations, we get. 


(2.55) 


f 


-1 

-1 

r 

‘‘■ch 


<N2>1 '«3>1 <N4>1 


(O' >- 

g 1 


1, = 

<Ni>2 (Nj)^ (N3)2 (N,)^ 

■ 

g 2 

’“-c's 


‘'^1*3 <«2’3 ‘'<3>3 <'^4>3 


g 3 



<Ni>4 (Nj), (N3)^ <N,)^ 


>A 

g 4 


« 


k a 


(2.56) 


This gives the stress component at 

all the 

nodes of 

the 

element. The 

state of stress at the corner 

of the 

neck 

IS evaluated 

by 

averaging the stress matrix 

over all 

the 

adjacent elements. 

Finally the 

equivalent 

stress 

i s 

calculated by 

the 

relation 




«j* = 

eq 

I 

2 

— + o» 


<2 

.57) 


I 

12 2 




where o'. and 


are the principal 

stresses. 

This is 

done 


to simplify the expression for the Von liises criterion. 


Now the criterion slates that yielding occurs when 
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^eq ~ ’ yiPld stress in tension. ( 2 . 58 ) 

Thus the numerical value of immediately tells u 

whether there it yielding or not. 
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CHAPTER III 

RESULTS AND DISCUSSION 

A computer program was developed to implement the 

analysis presented in chapter II, The program first 

calculates the stresses in the blade and the disc which 

provide the boundary conditions for the main problem. In 

the second stage, the main problem - analysis of the root 
region as a plane stress contact problem with friction 
IS solved using finite element method. The program takes 

the following data as input ! 

(1) blade, disc and contact geometry 

(2) material properties of the blade and the disc 

(3) rotational speed of the turbine. 

The output consists of 

(1) stresses and displacements fields at the contact 

interf ace 

(2) equivalent stress at the corner of the neck of the 

blade root. 

The data on the material properties of the blade 
and disc and the turbine speed used in the computations 
is as follows s 

5 2 

Young’s modulus, E = £ x 10 N/mm 

Poisson’s ratio, v - 0.25 

coefficient of friction 
between the blade and 
disc, ^®0. £0 



3fe 

of the blade 

and disc material, P = 8 x 10“^ Kg/mm^ 

speed of the rotor = 3 OOO rpm 

3.1 STRESSES IN TI-K: BLADE AND DISC : 

The stresses in the blade and the disc are 
calculated using the analysis of sections S.1 and 2.2 
respectively. The following geometrical data is used for 

these cal culalions . 


outer 

radius 

of the 

disc. 

^ - 

170 mm 

inne r 

radius 

of the 

disc. 

^2 = 

50 mm 

1 ength 

of the 

blade , 


'b = 

100 mm 


distance of the cross— section 
of the blade at which the 
axial stress is calculated, Z = 20 mm 
The axial stress on the blade cross-section, AB (Fig. 3.1) 

comes out to be 14.52 N/mm^. The variation of radial 

and circumferential stresses in the disc with respect to 

the radial co-ordinate is shown in Fig. 3.2. The stress 

component t^ and t^ on the boundaries CD, DE and EF 

(refer to Fig. 3.1) are calculated from this information. 
3.2 ANALYSIS OF THE ROOT REOICHnI : 

Analysis of the root region is carried out using 
the finite element formulation of section 2.3. Four noded 
isoparametric element is used. The program was tested by 


analysing 

the 

tee-root 

considered 

by Shrivastava 

(15) and 

comparing 

our 

results 

with 

his 

results. The 

comparison 

was made 

for 

the stresses 

and 

displacements 

only along 


the contact length and there was a very good agreement. 


i 



For the blade 
B = Width of the k 
L = Length of the shoulder 
T = Thickness of the shoulder 
£x= Clearance in x direction 
€z= Clearonce in z direction 

For disc 

C = Thickness of disc arm 


lL 


■TX 



AB and CDEF boundaries define the root domoin. 
Inner radius of disc = 50 mm. 


The X-root joint between the blade and disc 


CFrcumferentJal Stress MPo 
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Since no analytical or experimental results are available 
for the equivalent stress at the corner point, no such 
comparison could be made. 

3 . 2.1 RESULTS FOR ThC TEE-ROOT 

The following geometric data was used in the 
analysis of the tee-root (see Fig. 3.1) : 

B = IS mm 
C = 8 mm 
D = 5 mm 
L = a mm 


« =1 mm 

X 

tD = 1 mm . 
z 

The finite element mesh used in the analysis is 
shown in Fig. 3.3. For the tee— root the euivalent stress 

O' was calculated at the corner point A' (shown in Fig. 
eq 

3.1) for the values of T/L varying from 0.75 to 1.75. 


This was done keeping the value of L constant. The 

variation of o' with T/L is shown in Fig. 3.4. The 

eq 

graph shows that the equivalent stress decreases with the 


increase 

in 

the 

value of 

T/L 

follows . 

One 

can 

consider 

the 


beam subjected to distributed 
increased keeping L constant 
shoulders is increased. This 
Although the distribution of 
with T/L, total load remains 
change in bending moment at 
prominent as the change in 


. This can be explained as 
shoulder as a cantilever 
transverse load- When T/L is 
it means thickness of the 

increases its section modulus, 
the transverse load may vary 

the same. Therefore the 

the . corner A' may not be as 
section modulus with T/L. 


Since the section 


modulus increases with T/L, the 


bending 
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Fig. 3.3 Finite Element Mesh for Tee -Root Joint. 
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stress would decrease. Although & is not equal to 

e q ^ 

the bending stress, it is quite possible that a increases 

eq 

whenever bending stress increases. Thus the trend of 

Fig. 3.4 seems quite reasonable. 

3 . 2.2 RESULTS FOR THE STRADDLE ROOT ; 

Fig. 3.5 shows the geometry of the straddle root 
joint. The relations of various geometric parameters of 

the straddle root with the tee root are defined by 

eq.(1.1). These relations are reproduced here for the sake 
of convenience. 



As stated 

earlier. 

two different 

models 

of the 

straddle 

root are 

considered. 

In the first 

model 

thicknesses of 

both the 

necks are 

equal i.e. b^ 

= bg. 

In the 

second 

model the 

thickness 

of the lower 

neck 

is related 

to that 

of the first by the relation 
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Finite element mesh used in this analysis is shown in 
Fig. ( 3 . 6>a ) and (3.<bb). 

For the straddle root the variation of the 

equivalent stress at the corner points A' and B' (shown in 

Fig. 3.5) with respect to the ratio l^/l^ of the shoulder 

1 S 


lengths 

was 

determined. 

This 

was 

done for 

various 

values 

of T/L. 

As 

l,/l2 

increases , 

the 

length as 

well as 

the 

thi c)-.ness 

of 

the 

upper 

shoulder 

increases 

while for 

the 


lower shoulder these dimensions decrease. If these shoulders 
are considered as cantilever beam subjected to transverse 


distributed 

loads , 

then 

increase 

of 

l^/lg results in 

the 

increase 

of 

the 

section 

modulus 

of 

the upper shoulder 

but 

decrease 

in 

that 

of lower shoulder. 

It is possible 

that 


the bending moments at the corner A' and B' change with 

respect to However, the effect of the section 

modulus seems more dominant. Thus the bending stress and 

along with it the equivalent stress seems to decrease for 

the case of upper shoulder (point A' ) . For the lower 

shoulder (point B' ) it seems to increase. This trend is 

depicted in Figures (3.7) to (3.11). These graphs show the 

variation of the equivalent stresses at the two corners 


with respect to 

the 

ratio of 

shoulder 

lengths 

for 

five 

different values 

of 

T/L. The 

(a) part 

of 

the 

figure 


corresponds for the first model (b^ = 


bg) while the 


(b) 



AB onc3 CDEf boundones define the blode and disc straddle 


b ■= Width of the veck 
I = length of the sholder 
t = Thickness of the shoulder 
£x= Clearence in x direction 
£? = Clearance in z direction 


Subscript I to the upper 
shoulder 
Subscript 2 to the lower 
shoulder 



Fig 3.5 Straddle root joint between the blade 
and disc 




Fig.3.6 (b) Finite Element Mesh for Straddle root contact 
( Model n b| /b 2 = L / Ig ) 




part corresponds to the second model (fa^/bg = 1^/lg). 

The Figures 3.7 to 3.11 indicate the optimum value 

of which IS denoted by ^l-t/lg^opt' value 

at wfiich the equivalent stresses at the corner points 

become equal. For the given value of T/L, (1 /I ) is 

1 S opt 

lower for the second model. This means, in the optimum 
design, the second model has longer and thicker lower 
shoulder compared to the first model. This is probably 

required to compensate for the lower length of the lower 

necl of the second model. It may be recalled that the 

lower necl of the second triodel is thinner compared to the 
first model. The value of equivalent stress in the tee-root 
of the same T/L ratio is also indicated in the Figures 
(3.7) to (3.11). Figure (3.12) shows the variation of 

equivalent stress in the optimium straddle root with respect 

to T/L ratio. The variation shows the same trend as that 
of the tee-root. The values in the parenthesis indicate the 
value of ^^1^^2^opt‘ noted that, for given T/L 

ratio the second model has higher equivalent stress compared 
to the first model. 

3£.3 COMPARISION OF THE TEE AND OPTIMUM STRADDLE ROOT 

First we define the nondimensional equivalent 

stress in the optimum straddle root as follows s 


(O' ) . 

eq opt 


C(<y ) .3 

eq opt s 

(O’ )- 

eq T 


where (o' ) ^ = non-dimensional equivalent stress in the 

eq opt 

optimum straddle root 



96.D -1 


O' H 

94,0 - 
4 


o» the corresponding tee root is 106.2 MPc 



° Corner A 
* Comer 8 


110 - 

108 ' 

106' 

0 104 

1 

O-102 




T/L - 0.75 

of the corresponding tee root ie 

106.2 

MPa 


/* 

Q 

Comer k 



Comer B 

1 




Fig. 3.7 Variation of oquivolont stroseos at two condors 
with respect to the shoulder lengths 





Rq3.8 Variation of the equivalent stress at the two comers 
with re$p6ct to tho rotio of tho shouidftr lohgths*^ 



C4'- 




•j/lz. 

(b) b, / bj “ L/I 2 

Rq, 3.9 Voriation of the equivalent stress at the two comers 
wi^ respect to the ratio of the shouider lengthe-^ 


74^ T/L » 1.50 


of tne corresponding tee root is 72.05 MPo 



*1 /Ij 

(d) b,- b2 




eq MPo 


54 - 



T/L » 1.75 


gg T Of the corresponding tee root is 63.16 MPa 



(b) b, /b^« l/k 


Rq. 3.11 Variation of the equivoient stress at the two corners 
^ w'rth respect to the ratio of the shoulder lengths- 
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~ equivalent stress in the optimum straddle 
root 


^^eq^T ~ equivalent stress in the tee-root having 



the 

same T/L 

ratio. 




Variation of 

(& ) 
eq 

. with 
opt 

respect 

to T/L ratio is 

shown 

in the Fig. 

(3. 13) 

. The 

(a) 

part of 

the 

figure 

corresponds 

to the 

first 

model 

while the 

(b) 

part 

CO r res ponds 

to the 

second 

model . 

The figures 

are 

divided 


into two regions. In the region I, ^^eq^opt less than 

1. Therfore for T/L ratios lying in this region, it is 

advantageous to use the optimum straddle root rather than 
tee root. This is because it is subjected to lower 

equivalent stress compared to the tee-root of the same 
T/L ratio. In the region II, ^*^eq^opt greater than 1. 

This means for T/L ratios lying in this region, the 

tee-root has less value of equivalent stress compared to 

the optimum straddle root of the same T/L ratio. 
Therefore in this region, it is more advantageous to use 
the tee-root. It should be kept in mind that this 
comparison is only from the strength point of view and 

that too only for steady conditions. 

3.a4 CONTACT STRESSES AND DISPLACEMENTS : 

To show a typical variation of the contact 

stresses and displacements, the value of T/L 0.75 is 

chosen. At this value, the non-dimensional equivalent 
stress in the optimum straddle root is minimum. Figures 
3.14 to 3.17 show respectively the variation of normal 



5 7 


1 4 



T/L 


(a) b,» 



Bg.3.13 


on of non dimensionol «q 
:he optimum straddle root 


ulvalent stress 
wi^ respect to T/L 
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displace 

merit , 

tangential 

displacement , 

normal stress and 

the 

tangent 1 

a 1 s 

tress along 

the 

contact 


length for 

the 

tee-root 

. Fi 

gures 3. 1S 

to 3.21 

shows 

the 

variation of 

the 

same quantit 

les for the 

first 

model 


of the optimum 

straddle 

r 00 

t wh lie the 

Figures 

3.22 

to 

3.25 are for 

the 

second 

mode 1 

. In the 1 

ast two 

sets 

of 

figures, part 

(a) 

denotes 

the 

variation a 

long the 

upper 

sh 

oulder while 

part 

( b) IS 

for 

the lower 

shoulder. 

In a 

11 

three figures. 

X 

denotes 

the 

distance of 

the contact 

point from the 

axis 


of symmetry <Z-aj<is). These figures show the variation only 
for the right hand side contact surface, the normal 


stress 

and 

displacements 

are symmetric 

while the tangential 

stress 

and 

displacements 

are antisymmetric. 



The 

trend for 

the tee-root as 

well as both 

the 

mode 1 s 

of 

the optimum 

straddle root are similar. It 

is 

observed 

that the normal displacement 

decreases as we 

go 

away from 

the Z-axis. 

This indicates 

that the pull 

of 


the blade is more dominant nearer the Z-axis. This seems 
reasonable. Similarly the normal stress also decreases as 
we move away from the Z-axis. This can be explained as 
follows. For T/L = 0.75, thickness of the disc arm is 
more than the thickness of the shoulder (see Fig. 3.1). 
This makes the disc arm more stiff compared to the 


shoulder . 

Therefore the 

free end 

of 

the disc 

arm while 

deflecting 

applies more 

force on 

the 

shoulder 

on the disc 


arm. Tangential displacement increases as we go away from 
the Z-axis. The results on the tangential stresses don't 
appear to be very reliable although they roughly seem to 
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Rg 3.17 Variation of the Tangential Stress 

along the contact length (for Tee root) 
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Vcriotion of 
length 


(b) Lower shoulder 

tonqential displacement aion^g the contact 
(for straddle root, model D 



X co-ordinate mm 


(b) Lower shoulder 

FiQ.3.20 Variation of normo! stress along 
length (for straddle root. 


the contact 
model I) 
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Fiq.3.21 Variation of tangential stress along the contact 
length (for straddle root, model 1) 
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Rg.3.22 Variation of normal displocemsrrt olong the contoct 
length (for stroddle root, model II) 
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Rq. 3.24 Variation of normal stresses aiona th*. 

length (for straddle root, model in 
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pjn 95 Voriation of sheer stresses along the contact 
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1 nd 1 cate 

parabol i c 

variation. 

The 

erratic behaviour 

of 

the 

tangential 

stress 

may be due 

to 

the application 

of 

the 


slip boundary condition. This boundary condition is applied 
to the nodal forces rather than to the stress components 

on the contact surface. If it were possible to apply the 

slip boundary condition directly to the stresses, then 
probably the behaviour of the tangential stress would be 
more smooth. 
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CHAPTER IV 

CONCLUSIONS AND SUGGESTIONS FOR THE FUTURE WORK 

The following conclusions can be drawn on the basis of 

the results presented in Chapter III : 

(1) For the given T/L ratio, optimum straddle root is 

better if T/L ratio is roughly less than 1.5. This, 

however is only from the strength point of view and that 

too under steady conditions. When T/L ratio is roughly 

more than 1.5 tee-root should be used. 

(2) Among the models of the optimum straddle roots, the 

second model is slightly better as it has a lower 
value of equivalent stress. 

(3) Among various straddle roots of the given T/L ratio, 

the optimum straddle root is the most desirable as it 

has the same equivalent stress at the corner points. 

As mentioned in Chapter I, this work has several 

limitations in the form of idealisations made. In future, 
some or all of these can be accounted for, in particular 

the following modifications can be undertaken s 
(1) In case of straddle root, L should be varied to 

find the optimum ratio of shoulder lengths i.e. ^l-j/^g^opt* 
<2) The effect of nozzle forces should be considered. The 
results of this study may change completely with the 

inclusion of unsteady stresses. 

(3) Stagger of blade and radial placing of the root 
should be considered. Hence the root region should be 
analysed as a three dimensional contact problem. 
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